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Abstract 

We extend the framework of QCD factorization to non-leptonic B decays into three light 
mesons, taking as an example the decay B + tv + tv + tv~. We discuss the factorization 
properties of this decay in different regions of phase space. We argue that, in the limit of 
very large 6-quark mass, the central region of the Dalitz plot can be described in terms 
of the B —>► tv form factor and the B and n light-cone distribution amplitudes. The 
edges of the Dalitz plot, on the other hand, require different non-perturbative input: the 
B —)► 7V7V form factor and the two-pion distribution amplitude. We present the set-up 
for both regions to leading order in both a s and Aqcd/^ 6 and discuss how well the two 
descriptions merge. We argue that for realistic 5-meson masses there is no perturbative 
center in the Dalitz plot, but that a systematic description might be possible in the context 
of two-pion states. As an example, we estimate the B —)> pn branching fraction beyond 
the quasi-particle approximation. We also discuss the prospects for studies of three-body 
and quasi-two-body non-leptonic B decays from QCD. 



1 Introduction 


Three-body non-leptonic decays of heavy mesons constitute a large portion of the branching 
fraction. For B mesons, three-body non-leptonic branching ratios and CP asymmetries have 
been measured for a large number of channels, most notably by BaBar, Belle and LHCb |l|[ 6 j, 
and more is expected to come from LHC run 2 and from Belle II |7|. On the theory side, three- 
body non-leptonic B decays are interesting for several phenomenological applications, such as 
the study of CP violation and the extraction of the CKM angles a and 7 (see e.g. |8,9j). While 
in most cases there is a dominance of quasi two-body final states, in some decay channels 
the contributions from non-resonant three-particle states seems to be rather large |l0|. The 
study of the interference pattern of the resonances in Dalitz plots is a well established method 
to determine CP asymmetries 11 , 12 , while further information can be inferred on strong 


resonances, such as masses, widths and quantum numbers 13 


There are, however, two obvious problems in the quasi-two-body interpretation of resonant 
effects in multi-body decays, one practical and one conceptual. From the practical point of view, 
any parametrization of resonant structures is model dependent, as no universal line-shape for 
strong resonances is accurate, especially for broad states. On the conceptual level, the mere 
separation of resonant and non-resonant contributions is not clear-cut, most prominently in the 
case of non-leptonic decays where non-factorizable effects exist. 

In the case of two-body non-leptonic B decays, the heavy-quark limit has been exploited 
systematically in the context of QCD factorization 14-18 or Soft-Collinear Effective Theory 


(SCET) [19-23 , where the matrix elements factorize into a convolution of perturbative hard 
kernels, form factors and meson distribution amplitudes on the light cone. Corrections to fac¬ 
torization arise at subleading orders in the heavy-quark / large-energy expansion, and remain a 
source of uncertainty which is difficult to estimate. Some potentially important non-factorizable 
effects might be related to nearly on-shell intermediate states, such as charm-loops or rescat¬ 
tering effects from light mesons. Phenomenological investigations of such effects have limited 
potential, mainly because the kinematics of two-body decays is fixed. On the other hand, three- 
body decays have at their disposal a wide phase space where the energy dependence of such 
effects can be studied, with the potential of providing a deeper understanding of factorization 
and hadronic effects in B decays. 

It is fair to say that the theoretical description of three-body B decays is still in the stage 
of modeling. Common methods reflecting the state of the art are the isobar model (24j 25) 
and the K-matrix formalism j 261. In these approaches, resonances are modeled and the non¬ 
resonant contributions are often described by an empirical distribution in order to reproduce 
the full range of the phase space |27|. In the context of factorization, in Refs. 28-30 the ma¬ 


trix elements were factorized naively and the resulting local correlators were computed in the 
framework of Heavy-Meson Chiral Perturbation Theory (HMChPT), but no attempt was made 
to address the breakdown of factorization or HMChPT in the respective regions of phase space 
where they are not expected to apply. Other recent work relying on pQCD |3l], seems to repro¬ 


duce experimental values for CP asymmetries integrated in certain regions of phase space 32 


However, if the conceptual issues regarding the pQCD approach 33,34 cannot be resolved, 
its predictive power remains limited. In the future, novel model-independent approaches that 
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directly access CP violation (such as the Miranda procedure 112,35)) or methods based on fla¬ 


vor symmetries 36,37| could become interesting; however, for a quantitative description of the 


differential Dalitz distributions including amplitude phase information, a QCD-based approach 
is unavoidable. 

In the present letter we take a step in this direction, and study the factorization properties 
of charmless three-body and the corresponding quasi-two body 5 decay fQ For that purpose, we 
focus specifically on the decay B + —? 7r + 7r - 7r + , assuming that the 6-quark mass is large enough. 
We start by identifying the different regions in the Dalitz plot where the well-established fac¬ 
torization properties of two-body decays apply to the three-body case. In the heavy-mass limit, 
we discuss how to compute the central region of the Dalitz plot as well as its edges. We will 
see that the methods and the theoretical inputs are different in the different regions: while the 
center can be described in terms of regular from factors and pion distributions, the descrip¬ 
tion at the edges requires introducing generalized versions of these hadronic matrix elements. 
Generalized form factors and distribution amplitudes have been already studied, the former in 
the context of semileptonic B decays |40|, and the latter in connection with two-meson electro¬ 
production |41,42 or semileptonic r decays (43) (see e.g. 44-51,53,54,63)) As an application, 


we consider the B + —> p° 7r + branching ratio by integrating the differential rate around the p 
resonance. Finally, we discuss how both descriptions merge to describe the full Dalitz plot, and 
what we can expect for realistic 6-quark masses. 


2 Identifying regions in the Dalitz plot 

We consider the decay B + —? 7r + 7r - 7r + , and define the external momenta as: 

B + {p) —> 7r + (fci) 7r _ (£; 2 ) n + (kz) with p = k\ + fc 2 + an d E± M < E^ M , (2-1) 

where CM refers to the 5-meson rest frame. We neglect pion masses in the kinematics, such 
that: 

p 2 = m 2 B , k 2 = 0 , Sij = — + 2 ^ = — 2 fcj (i ± j ) . (2.2) 


rriz 


l b m B 

The kinematics of the three-body decay is completely determined by two of the three kine¬ 
matic invariants Si 2 = Si 3 = s ++ and s 23 = .s l |A h , which (in the massless limit) satisfy 
S 12 + S 13 + S 23 = 1 and 0 < Sij < 1. The physical kinematical region in the plane of two 
invariants (the Dalitz plot) is given in this case by a triangle (see Fig. [I]). 

We distinguish three special kinematical configurations: 

I. “Mercedes Star” configuration: This corresponds to the central region of the Dalitz 
plot, where all the invariant masses are roughly the same and of order of mg: 


Region I : s ++ ~ s 1 ™ ~ ~ 1/3 


(2.3) 


corresponding to the kinematical situation where all three pions have a large energy in 
the 5-meson rest frame and none of the pions moves collinearly to any other. 


1 The main ideas developed here have been discussed qualitatively by M. Beneke 38 and I. Stewart 39 
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Figure 1: Left: The physical kinematical region in the plane of two independent momentum 
invariants s 1 ”™, (Dalitz plot), divided into the different regions with special kinematical 
configurations: I - Mercedes Star configuration, IIa,IIb - Two collinear pions, IIIa,b - One soft 
pion. Right: Dalitz plot distribution for B + —> 7r + 7r - 7r + from Ref. |5j 


II. Collinear decay products: This corresponds to regions of the Dalitz plot where one 
invariant mass is small and the other two are large. The kinematic configuration is such 
that two pions are collinear, generating a small invariant mass recoiling against the third 
pion. In our case there are two such regions: 

Region Ila : ~ 0 , s 1 ^ ~ ~ 1/2 (2.4) 

which is the region where the two 7r + move collinearly, recoiling against the ir~ , and 

Region lib : s 1 ™ ~ 0, s ++ ~ s+£ h ~ 1/2 (2.5) 

where the ir~ and one 7r + move collinearly, recoiling against the second 7r + . 

III. One soft decay product: The regions of the Dalitz plot where two invariant masses 
are small and one is large correspond to kinematical configurations where one pion is soft 
and the other two are fast and back-to-back. In our case there are two such regions: 

Region Ilia : s ++ ~ s 1 ™ ~ 0, s+£ h ~ 1 (2.6) 

which is the region where one 7r + is soft, and 

Region Illb : s 1 / 1 ? ~ .s 1 /^ ~ 0, s ++ ~ 1 (2.7) 


where the 7r is soft. 

The different regions are shown in Fig. [l} For a very heavy B meson, region I is dominant, 
since the condition m 2 B Sij Aq CD is satisfied in most of the Dalitz plot. The edges of the 
Dalitz plot, corresponding to collinear and soft configurations, will be small. However, in the 
edges, all the resonances show up, corresponding to quasi-two particle decays. The masses m,R 
of these resonances do not scale with the heavy b quark mass and hence the width of regions 
II and III scale as mR/mB , showing the dominance of region I in the infinite mass limit. 
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Figure 2: Factorization formula in the center (region I) 

In the following we propose to perform a QCD factorization calculation in region I of the 
Dalitz plot in terms of the pion and 5-meson light-cone distributions and the B —>■ n form 
factor. The presence of resonances in region lib will be signaled by a singular behavior of 
the factorized amplitudes of the form 1 /.s 1 ° w \ A proper treatment this region requires to set 
up a different calculational method, which requires a different form of QCD factorization. In 
this case, new non-perturbative quantities need to be defined: a light-cone distribution for two 
pions, and a form factor for the B —>■ mv transition. Finally, it is important to check that the 
two calculations match properly in order to obtain a complete description of all regions in the 
Dalitz plot. 

The region Ha contains no resonances (corresponding to the 7r + 7r + channel), and will see 
that the factorized amplitudes are regular as s ++ —>■ 0. In this case the “perturbative” result 
should provide a good description of the rate when integrated (or smeared) over a suitable 
interval, in the sense of parton-hadron duality. 


3 The central region of the Dalitz plot 


Our starting point is the heavy-quark limit, where we assume that m^/s /3 Aqcd- In the 
central region of the Dalitz plot (region I) we have all invariant masses of the order 3 and 

hence we expect the factorization formula 


(7T 7T 7T 


-\Oi\B+) s 


-1/3 


= T 


?B—>77 


4, ® 4„ + T‘‘ ® 4 B ® 4„ ® 4„ ® 4, 


(3.1) 


where Oi is a four quark operator in the effective weak Hamiltonian. This factorization formula 
is illustrated in Fig. [2} The hard kernels can be computed perturbatively in QCD. Some typical 
diagrammatic contributions are shown in Fig. [3} We will consider here only the leading a s 
corrections, and neglect next-to-leading a 2 S contributions such as (c) and (d) in Fig. |3j While 
the study of c? s corrections is beyond the scope of this analysis, we expect these to be about 
~ 10% relative to the leading color-allowed amplitude, similar to the case of B —» 7T7t (see 
e.g. Ref. 155]). The diagram (b), where the gluon is ejected from the spectator, requires the 
spectator quark in the B meson to have a large virtuality of order rn h , which is either suppressed 
in the heavy-quark limit or requires an additional hard interaction. All in all, we do not include 
the second term in Eq. (3.1) in our analysis, nor radiative corrections to Tj. To this order, 
the convolutions of the hard kernel T/ with the B —> 7r form factor and the pion light-cone 
distribution can be computed without encountering end-point singularities. While this would 
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Figure 3: Sample contributions to the hard kernels Zf and Tj 1 in the factorization formula for 
the central region. The leading a s contributions are given by diagrams such as (a) and (b), 
while (c) and (d) are next-to-leading in a s . See the text for details. 


be a trivial statement in the case of two-body decays, we stress that here the kernels T*(u,v) 
already depend on the momentum fraction of the quarks at the leading order, making the 
convolutions non-trivial. 

The differential decay rate d 2 T/(ds ++ ds + _) computed in this way shows some interesting 
features. First of all, moving from the central point s ++ = s + _ = 1/3 (region I) toward the 
edge s ++ 0 (region Ha), we find that the rate remains regular, that is, it approaches a finite 

limit as s ++ —y 0. This can be seen explicitly in the calculation, with no propagator becoming 
soft as —>■ 0. More precisely, moving away from the center along the line s + _ = (1 — s ++ )/2, 

we find 


dr 

ds .|_ + ds + _ 


S++ ■ 

S+- 


► 0 
1/2 


r» /+(m|/2 ) 2 


(3.2) 


up to a coefficient of order one, with 

r "“ 32^ ' 

Here f+(q 2 ) denotes the vector B —>■ n form factor, defined as: 


(3.3) 


q,u «P ~ q)\b*q\B(p )) = (m| - ml)f 0 (q 2 ) ~ (m| - q 2 )f + (q 2 ) , (3.4) 

where in the last term we have employed the large-recoil-energy relation (56 . 

In Fig. |4] (left panel) we show the exact dependence of the rate as a function of s ++ , along 
this direction in the Dalitz plane. This regular behavior does not depend on how we approach 
the s ++ = 0 edge. 

The situation is very different if we consider the behavior of d 2 T/(ds ++ ds + _) as s + _ gets 
small (towards region lib). We consider now the direction along the line s ++ = (1 — s + _)/2. 
In this region the rate behaves as, 


dr 

ds .|_ ds .|_ 


s -\— —> 0 
S++ ~ 1/2 



U{m 2 B /2f 


regular terms as s + _ —>■ 0 


(3.5) 


rendering the rate non-integrable. This behavior is expected, as the edge of the Dalitz plot with 
small s_|_is determined by hadronic resonances, dominantly the p resonance. In this region 
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Figure 4: Differential decay rate when extrapolated from the center region of the Dalitz plot 
towards the collinear edges. The extrapolation to small s ++ remains regular, while the limit to 
small s + _ diverges. See text for details. 

the three-body decay effectively becomes a quasi two-body decay, and the methods to describe 
collinear and soft parts of the Dalitz plot have to be modified. 

4 The collinear regions of the Dalitz plot 

The case where one invariant mass is small is kinematically very similar to a two-body decay. 
We expect then a similar factorization theorem, with the difference that one of the particles in 
the two-body case is substituted by a pair of particles with small invariant mass, which must 
be described collectively. In order to describe the soft and collinear regions of the Dalitz plot, 
we therefore need to introduce additional non-perturbative quantities: the two-pion light-cone 
distribution (27 tLCD) amplitude and the B —» tttt from factor. 

To leading twist, the 27 tLCD for a 2-pion system (7r + 7r - ) is formally given by the matrix 
element [41 

Slp(z, h,k 2 ) = hjf J dx~e~ lz ^ x ~ )/2 ('K + (k 1 )TT~(k 2 )\q 0 (x)[x,O]q a (O)\O) x+=x±=o (4.1) 

where {a, /?} are Dirac indices, q = u,d, and [x, 0] is a Wilson line. We take k l2 = k\ + k 2 and 
define two light-like vectors n± = (1, 0, 0, ±1) such that 

k u = ^ n + + ^ n - and x M = y u + + Y n ~ + ’ ( 4 ' 2 ) 

and such that when k\ 2 —> 0, then ky 2 —> 0. The variable z is the fraction of the momentum 
h '12 carried by the quark q. 


42,46,48 
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The Lorentz decomposition of the matrices S'‘ yJi consistent with parity invariance, keeping 
only terms that contribute at twist-2, is given bj|^] 

S q a p = ^$J(^,C^)^12 + $1(aC,^)W^^2 , (4.3) 

which defines the vector ($||) and tensor (<E>_i_) 27 tLCDs. The variable C = /^i 2 is the light- 

cone momentum fraction of 7r + . In terms of invariants, we have 

k\ 2 = m 2 3 B s \2 , (= - — . (4.4) 

J- — S 12 

Isosinglet (d>° = 4[$ u + $ d ]) and isovector (4> x = 4[$ u — 4> d ]) 27 tLCDs have been discussed 
in the literature (e.g. Refs |46][48j). The vector 1=1 27tLCDs are normalized as 

J dz<S>l(z,(,s) = (2C-1 )F«(s), (4.5) 


46 


where F n (s) denotes the pion vector form factor. C-parity and isospin invariance imply that the 
corresponding integral is zero for the isosinglet component |46|. At the leading order, the hard 
kernel T(z) with which (, s ) is convoluted in the B —> 7r7T7r amplitude does not depend on 

the momentum fraction z, so the amplitude depends only on the local form factor F n (s), just 
as the leading contribution in B —» 7T7t depends only on f n . In addition, at the leading order 
the tensor distribution 4 >t does not contribute. 

The vector form factor F 7T (s) in the time-like region (s > 0) can be obtained from measure¬ 
ments of the process e + e _ —> |57(^]- see Fig. [5| We employ here the fit parametrization 


2 This definition of agrees with reference 46 (up to isospin decomposition, see later). However the 
definition for <l>x might differ from that in 46 by an overall factor, which we do not address here because to 
the order considered 4>x will not appear in the amplitude. 


3 For simplicity we use only the latest Babar data, but see also Refs. 


58 


64 
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of Ref. 165], which is consistent with general principles of QCD at low energies, and covers the 
energy range of interest, including the relevant resonances in that range. The particular choice 
of parametrization is not very important for the absolute value |i^(s)|, where a good fit to the 
data is enough (see Fig. [5]), but it is important for the phase, where data is not so precise. A 
thorough analysis of the phase of F n (s) and its impact in B —> mm is beyond the scope of this 
paper, but it becomes a crucial issue as soon as one attempts to describe CP asymmetries. We 
leave this for future work. 


The second nonperturbative input is the B —» mi form factor, which has been discussed 
already in the context of B — > rnifv decays in Ref. |40|. We consider the generic form factor: 

F a p(k 1 ,k 2 ,k 3 ) = ('K + (k 1 )7r~(k 2 )\b /3 u a \B + (p)) , (4.6) 

where a, (3 are Dirac indices. The most general Lorentz decomposition consistent with parity 
invariance is given in terms of four independent form factors: 


F = F t 


r#375 + F 2 ^(o)75 + F 3 #(||)75 + F4 e a /3 7/i A;“ kj^ 


-F t 75 , (4.7) 


! 4yAf '4 (m b + m u y 

with fc(o)j^(||) two orthogonal space-like vectors. Due to the structure of the leading order 
contributions, the time-like form factor F t (C,s 12 ) will be the only one relevant here. This 


definition of F f coincides with Ref. 40 


In order to be able to make a quantitative prediction, we can relate the different B —> mr 
form factors to the 27 tLCDs via a light-cone sum rule |66|. For the time-like form factor F t we 
liavcEl 


m,su) = 


mi 




dz 

— exp 

z 


(1 + s 12 z)m% rn 


,2 1 


M 2 


zM 2 




(4.8) 


no 


where fs is the static R-mcson decay constant extracted from a corresponding sum-rule, which 
is correlated to the Borel parameter M and to the threshold parameter u 0 . These three param¬ 
eters must be determined simultaneously with the condition that the physical decay constant 
and form factor are independent of M and uq. While we do not attempt to perform a full error 
analysis here, we note that the values }b — 0.316, Uq — 0.6 and M 2 ~ 10 GeV 2 satisfy this 
correlation approximately. In the asymptotic limit, given by <3>jj = 62(1 — z)(2( — l)F n (s) 
and setting \Jk f = m-. we have 

1 


48 


F t 1 (C,s 12 ) = 


3^2m 2 (2C-l)^( gl2 ) 
fBm w 


dz z exp 


(1 + si 2 z)m 2 B m\ 


,2 1 


M 2 


zM 2 


(4.9) 


no 


With this we have all ingredients for the factorization formula valid in the collinear regions 
of the Dalitz plot. The modified QCD factorization formula reads, in terms of the new non¬ 
perturbative quantities: 


(^ a ^ b ^ c pi\B) Sab ^i = 

+ T n 


?B—>7T C 


•F-a-t + T. 


I 

ab 


jB^7T a 7T b 




® $71 


^ 7T a 7T b 


(4.10) 


This is a tentative expression where we have ignored a possible contribution from the distribution We 


use this formula for illustrative purposes. The final form of this expression will be presented in Ref. 66 





















la 

n 


7T C 


Figure 6: Factorization formula for the edges of the Dalitz plot (region II). 


This formula is illustrated in Fig. [6] and yields now the description of the Dalitz plot in the 
kinematic regions I la and lib in Fig. [TJ 

Using the QCD sum rule relation ( |4.8| ) for the B —> 7T7t form factor, and writing the 27 tLCD 
in terms of the pion form factor F nj we may now write down the amplitude in the region of 
small s_|_, at leading order and leading twist. We find Q 


*4|s+- 


S + -<1 


= —\f 

V2 [h 


m, 


-( a i — a i) • Ft(C, s+-) + m 2 B (a 2 + a4)(2^ — 1) • /o(s+-) • ^r( s +-)] > (4-11) 


where the parameters a; are combinations of Wilson coefficient^] (see e.g. Ref. fl6|). To this 
order all the convolution integrals are trivial. Again, we neglect a s corrections and hard- 
scattering with the spectator quark for simplicity. Hard kernels are known already at NNLO 


from studies of two-body decays [67- 73 , but the convolutions with two-pion distributions still 
need to be worked out. In particular new distributions appear (e.g. <I>_l) that do not contribute 
at the leading order. This is beyond the scope of this work. The conclusions derived here at 
this order of approximation should nevertheless remain valid. 

A qualitative difference of three-body decays in this kinematic regime with respect to two- 
body decays is that the nonperturbative input is much richer in terms of QCD effects. In 
particular, F n contains resonance and rescattering contributions, including an imaginary part 
from non-perturbative dynamics, in contrast to two-body decays where strong phases are, at 
the leading power, of perturbative origin. This has implications both for quasi-two-body decays 
and for CP asymmetries. Most of the information on F 7T (s) can be obtained from data (see 
Fig- iu, allowing for a data-driven model-independent interpretation of three-body Dalitz plots, 
at least within the accuracy of factorization theorems. 

As a simple application of this result, we estimate the branching fraction BR(B + —» pi r + ) 
by integrating the differential decay rate in a neighborhood of the p resonance: 


BR(B + —y pir + ) = / ds ++ / ds + - 


t b dl 
ds ++ ds + - 


= / ds 


++ 


ds+- 


T B m B \A\‘ 

32(2tt) 3 


(4.12) 


5 Here F t corresponds to the combination F t ° + . The convolution with <hj| in the sum-rule for F t ° is in 

general not zero because the integrand is not even in z. Since the isosinglet distribution is mostly unknown, even 
in the asymptotic limit, we will nevertheless disregard this term altogether in the numerical analysis, keeping 
in mind that this issue requires further investigation. 

6 More specifically, we have ai ? 2 = V* b V u d(Ci j 2 + C 2 ,i/N c ) and a 3,4 = V* h Vtd{Cz^ + C^z/N c ), with the Wilson 
coefficients C« defined as in Ref. 


16 
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Figure 7: Differential decay rate obtained from the description in terms of two-pion distributions 
for small s + _. Left: extrapolation to s + _ ~ 1/3, with the p"(1700) apparent, and the p — oo — p' 
peak, in logarithmic scale. Right: Zoom to resonant contribution from the p(770) and w(782). 


where s± = ( m p ± S) 2 /m 2 B , and we will take S = nT p , with n specifying the cuts in units of the 
p-meson width. We find: 


BR(B + -> pvr + ) ~ 2.4 • 10“ 6 

for 

n = 1 

(4.13) 

BR(B + ->■ p7T + ) ~ 3.0 • 10“ 6 

for 

n = 2 

(4,14) 

BR(B + ->• pvr + ) ~ 3.2 • 10“ 6 

for 

n = 3 

(4.15) 

BR(B + -> pvr + ) ~ 3.3 • 10“ 6 

for 

n = 4 

(4.16) 


We note that extending the cuts beyond m p ± 4r p does not modify the result very much, as 
the resonant p contribution dominates the full decay rate. Comparing these numbers to the 


experimental value 74 


BR(B + ->■ p7r+) exp = (8.3 ± 1.2) • 1(T 6 (4.17) 

we see that the result is in the right ballpark. However, BR is an object different from the 
B —>■ p7r branching fraction as given in |T4], and can be measured experimentally in a direct 
and model-independent manner, without the need to extract the p from the full distribution. 
At this point we must emphasize that this is still a very crude estimate, and a more careful 
study would need to be performed to really test the data. 


5 Discussion 

So far we have used two different factorization formulas for region I and region II. Region I 
has been described using the conventional QCD factorization in terms of single pion states 
(which we will call QCDFj hereafter), while region II has been described in terms of hadronic 
input describing two-pion states with small invariant mass (called QCDFn hereon). To get 
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the full Dalitz distribution one needs to match the result from the central region with the one 
of the edges. To this end, we assume that there is an intermediate region between the edge 
(s 1 ™ = s ~ 0) and the center (.s/A = s ~ s+!l h ~ 1/3) where both descriptions apply. This 
region corresponds to ti^ CD /m 2 B <s< 1/3, and it certainly exists if mg is large enough. We 
will investigate below whether this happens for realistic 5-meson masses. 

In this intermediate region, one might use QCDF n (as in Section |4j) to write the amplitude 
in terms of two-pion states, then take the perturbative limit for the 27 tLCDs and B —>• mr form 
factors, and finally compare the result with the factorized QCDFi amplitude of Section [3| The 
idea is that, for s Aq CD /m£, we have (schematically 

^ttt r -t fl J dudvT^u.v)^^)^^) , (5.1) 

F b ^ f n F B ^ n (0) J duT F (u,v) cf) n («)+••• . (5.2) 


Taking this limit for the leading power contribution in QCDFn, one recovers fully some of the 
contributions obtained using QCDFi. 

In Table [l] we show the correspondence between the different contributions to the ampli¬ 
tude in this intermediate region, either in QCDFi or QCDFn- The first column shows the 
contributions from two-pion distribution amplitudes. In QCDFn (lower diagram), this is a 
leading-power contribution proportional to the 27 tLCD, $ W7r . As the invariant mass of the 
two pions in this intermediate region is also large, the two pions can be factorized according to 
Eq. (5.1). The production of two pions with large invariant mass requires a hard gluon, as shown 
in by the diagram at the top (corresponding to QCDFi). A similar argument goes through for 
the B —>■ 7T7T contribution, shown in the second column. The contribution in QCDFi (where 
the two pions are assumed to have large invariant mass) requires a hard gluon (top diagram), 
and can be obtained from the contribution in QCDFn (bottom diagram) by factorization of 
accorc ii n g to Eq. (5.2). We have checked analytically that applying Eqs. (5.1), (5.2) to 
the amplitude in Eq. (4.11) we recover the corresponding results in Section [3j 

However, some contributions in QCDFi correspond to contributions in QCDFn that are 
power suppressed, and do not arise from the perturbative limit of leading power contributions 
in QCDFn- These are shown in the last four columns in Table [lj Again, the contributions in 
QCDFi (s Aqcd /%) require a hard gluon. Columns 3 and 4 show the cases in which this 
gluon becomes collinear (in the [7T7 t] direction) as s —>• 0. They are termed “non-factorizable” 
since the gluon connects the two different collinear sectors. As s —>■ 0, the quark propagator 
remains hard, which represents a power suppression with respect to the leading contributions. 
Columns 5 and 6 show the cases in which the gluon remains hard for all s < 1/3. For s —> 0, 
these match onto 6-quark operators that are again power-suppressed with respect to the leading 
contributions. There is therefore a one-to-one diagrammatic correspondence between QCDFi 
and QCDFn, but this correspondence does not respect the power counting. 

We note at this point that in the center, since all invariant masses are large and of order 
, there are always two hard propagators, leading to an amplitude that is power suppressed 


m 


B 


7 The factorization of 2-7 tLCDs in the perturbative limit has been studied in Ref. 
B —> 7T7T form factors will be discussed in 76 


75 


The factorization of 


The dots in Eq. (5.2) account for “factorizable” contributions 


proportional to the R-meson light-cone distribution, corresponding to neglected contributions in Section [3] 
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Table 1: Diagrammatic correspondence between the different contributions in QCDFj and 
QCDFu. Crosses denote alternative insertions of the gluon. One- and two-pion distributions 
are denoted by (j) n and <p nn respectively, while F n and F n7r denote B —>■ it and B — > tttt form 
factors. The last four contributions are leading at the center but power-suppressed at the edge. 
See the text for details. 


with respect to the amplitude at the edge. In addition, the perturbative nature of the hard 
gluon exchange leads to an a s (rrib) suppression at the center, which is not present at the edge, 
where the gluon becomes soft. All in all, the amplitude at the center is expected to be both 
power- and aq-suppressed with respect to the amplitude at the edge. 

While the previous considerations imply that formally there must be a good matching 
between both regions, the question is whether this happens in practice for realistic 5-meson 
masses. To this end we focus on the 27 tLCD contribution shown in the first row in Table [H 
This contribution arises from the second term in Eq. (4.11). We find that, in the limit of large 
(m B s .(__), this amplitude reproduces the corresponding contribution obtained from the QCDFi 
calculation in Section |3|P| The particular values of s + _ for which this matching occurs depends 
on the value of m 2 B . In Fig. [8] we show the results of both calculations for different values of 
ms■ We see that for rn B ~ 20 GeV there is enough phase space to reach a perturbative regime 
in the central region of the Dalitz plot. However, the phase space gets reduced considerably 
when ms is decreased to its real value, where there seems to be no perturbative regime, i.e. 
the Dalitz plot is completely dominated by the edges. 

Similar conclusions are expected for the B —>• 7T7t form-factor contribution in the second 
column in Table [lj Adding the rest of the central-region contributions to the perturbative 
side, we will get a mismatch at large ( m 2 B s + _) of the order of the perturbative contribution 


8 This happens by construction, since we force the function F~(s) to satisfy the perturbative limit asymptot¬ 
ically for large sm 2 B . This is in fact the only information we have on F 7r (s) at large energies, since data reaches 
only up to ~ 3 GeV. For our purposes, the relevant observation is that data shows that the perturbative regime 
might lie beyond 3 GeV. 
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Figure 8: Contributions from 27 tLCDs to the B + —y differential branching fraction, 

for s ++ = (1 — s + _)/3: Full contribution (solid) and perturbative contribution (dashed). A 
perturbative region exists for large s + _ in the heavy-quark limit, but probably not for realistic 
values of the 6-quark mass. 

itself, which is expected to be of the same order as neglected power corrections to the QCDF n 
calculation in the perturbative limit. This corresponds to the contributions in the last four 
columns in Table [l] Since we anyway do not expect the two-pion system to factorize into 
single-pion distributions as early as m 2 B s + _ ~ 8 GeV 2 , we can conclude that the QCDFi 
calculation of Section [3] might not be relevant in any region of the Dalitz plot. 

6 Summary and conclusions 

Three-body B decays provide many opportunities for studies of flavor physics and CP viola¬ 
tion, as well for studies of factorization issues in QCD. While a large amount of experimental 
information is already available, an even larger amount is expected from future studies at the 
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LHC, as well as from Belle-II. These promising experimental prospects are not yet backed-up 
by theoretical studies able to describe these decays differentially in the kinematics, and in a 
model-independent manner. In this paper we have performed a first study in the context of 
QCD factorization. 

The Dalitz plot for three-body B decays can be divided in several regions with special 
kinematics and with different factorization properties. In the heavy-quark limit, the amplitude 
in the central region of the Dalitz plot factorizes into regular B —y tt form factors and pion light- 
cone distribution amplitudes. On the other hand, near the edge two pions become collinear 
and the amplitude resembles a two-body decay, with the difference that a new type of non- 
perturbative functions must be introduced: B —>■ nn form factors and 27 tLCDs. The fact 
that these objects cannot be factorized further is signaled by the divergence of the factorized 
expressions at the center when one invariant mass is taken to zero. We have calculated the 
amplitudes in both regions at the leading order, and verified these factorization properties. 

Assuming that these two regions are well described by the respective calculations, it is 
interesting to determine how well these two descriptions merge at intermediate kinematical 
regimes. We have seen that some of the contributions at the center correspond, in the heavy 
quark limit, to the expression for the amplitude at the edge with factorized B —» 7 T 7 t form factors 
and 27 tLCDs. Therefore, a parametrization of these nonperturbative objects that is consistent 
with their perturbative limit leads automatically to a well behaved limit of the result at the 
edge when extrapolated to the center. However, it seems that the perturbative regime is only 
kinematically allowed for 6 -quark masses several times larger than the real value. In addition, 
the rest of the QCD factorization contributions at the center correspond to power-suppressed 
corrections at the edge. 

At this point, a more refined study of the Dalitz plot distribution based on a factorization 
in terms of two-pion distributions and B —> 7nr form factors seems worthwhile. Next-to-leading 
(NLO) corrections to the hard kernels are already known from two-body decays, and can be 
used directly in the factorization formula of Eq. (4.10) to verify factorization at NLO. On the 
phenomenological level, this requires a better knowledge of two-pion distributions. First, one 
must go beyond the local limit (where information other than F n (s) is needed). Second, the 
tensor distribution is expected to contribute, while little is known about it. Besides the 
traditional studies of two-pion distributions from 7*7 —> mr and r —> mrv T1 we propose to 
study such distributions in the context of non-leptonic decays such as B° s —> D~ 7 r + 7 r° where 
factorization is considerably simpler than for B —» it nit. Regarding B —* 7 T 7 t form factors, 
better knowledge is also required. Improved light-cone sum-rule calculations [ 66 ], as well as 


precise experimental studies of the semileptonic decays B —y mriis 40 and B —y tttt/T will be 
essential. 

Besides B —» tttttt decays, other three-body decays with kaons (B + — > K + ir~ir + , etc.) have 
been studied experimentally at A-factories and the LHC (e.g. [|, 6 i|). Their branching fractions 
are higher because they are not CKM suppressed, with the corresponding impact in terms of 
statistics. These channels can be studied in a similar fashion. This requires knowledge on 
B —» Kn and B —>■ KK form factors, as well as Kn and KK distributions. Again, these can 
be accessed from semileptonic B decays (e.g. |77|) and r decays (e.g. |78|). 
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